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We present a lattice Boltzmann study of the hydrodynamics of a fully resolved squirmer, confined 
in a slab of fluid between two no-slip walls. We show that the coupling between hydrodynamics and 
short-range repulsive interactions between the swimmer and the surface can lead to hydrodynamic 
trapping of both pushers and pullers at the wall, and to hydrodynamic oscillations in the case 
of a pusher. We further show that a pusher moves significantly faster when close to a surface 
than in the bulk, whereas a puller undergoes a transition between fast motion and a dynamical 
standstill according to the range of the repulsive interaction. Our results critically require near-held 
hydrodynamics; they further suggest that it should be possible to control the density and speed of 
squirmers at a surface by tuning the range of steric and electrostatic swimmer-wall interactions. 

PACS numbers: 47.63.mf, 87.17.Jj, 47.63.Gd 


Introduction: Motile organisms such as bacteria and 
sperm cells have a natural tendency to be attracted to¬ 
wards surfaces, and to swim near them [1]. This phe¬ 
nomenon may be relevant for the initial stage of the for¬ 
mation of biofilms, the microbial aggregates which often 
form on surfaces. Experiments have shown that this ten¬ 
dency is not unique to living swimmers, and is also exhib¬ 
ited by phoretic, synthetic active particles laa]: in that 
context, it has been exploited for example, to attract mi¬ 
croswimmers inside a colloidal crystal, where they orbit 
around the colloids |1]. The interaction between self- 
propelled particles and walls also provides a microscopic 
basis for the rectification of bacterial motion by asym¬ 
metric geometries {e.g. funnels) [5]. 

Previous work has proposed two possible mechanisms 
for surface accumulation of self-motile particles. A first 
view is that accumulation occurs through far-held hy¬ 
drodynamic interactions [3]. Another possibility is that 
motility itself, in the absence of solvent-mediated inter¬ 
actions, leads to accumulation [ZlIH]: this mechanism re¬ 
quires a small enough channel, where the gap size is of the 
order of the typical distance travelled ballistically by the 
active particles, before rotational diffusion or tumbling 
reorients them. The case of phoretic particles may be 
more complex mu, and may depend on the dynamics of 
the chemicals reacting at the swimmers surface. Schaar et 
al. recently showed that hydrodynamic torques strongly 
affect the “detention times” over which microswimmers 
reside near no-slip walls m- The behaviour is partly 
controlled by the details of the force distribution with 
which active particles stir the surrounding fluid; previous 
work has also shown that these determine the equilibria 
of a spherical squirmer near a flat surface mma. 

Previous theories have not systematically studied the 
effects of a short range repulsion between the particle and 
the wall; in practice this is always present in experiments. 


either due to screened electrostatics, for charged walls, 
or due to steric interactions, e.g. for polymer-coated sur¬ 
faces. Here we show that explicitly including this re¬ 
pulsion is important, and strongly affects the dynamics 
near a surface. The interplay between hydrodynamics 
and short range repulsion can lead to trapping, periodic 
oscillations, and to a swimming speed significantly dif¬ 
ferent from that in the bulk. These results provide an 
experimentally viable route to tune microswimmer con¬ 
centration and speed near a no-slip surface. Furthermore, 
recent theoretical calculations m predict that the equa¬ 
tions of motion for squirmers which are pushers (exerting 
extensile forces on the fluid) do not possess any station¬ 
ary bound state solution, i.e., where the particle swims 
stably along the wall at fixed orientation. Our finding of 
oscillatory near-wall dynamics shows that even without 
such stationary solutions, trapping of swimmers at walls 
is possible. 

Squirmer model: A popular model for swimmer hydro¬ 
dynamics is a spherical squirmer - a particle which is 
rendered self-motile through a surface slip velocity m- 
The squirmer model has been used to study, e.g., the 
collective motion and nutrient uptake of swimmers in 
thin films [Ill US] ■ To study the dynamics of a model 
squirmer in confinement we employ a lattice Boltzmann 
(LB) method |16j . To achieve time independent squirm¬ 
ing motion, the following tangential (slip) velocity profile 
at the particle surface is used m 


u{e) = 2j2 

n—1 


sin0 dPn{cos9) ^ 

n{n + l) dcosO 


( 1 ) 


where 9 is the polar angle and is the nth Legendre 
polynomial mi. 

In the LB method a no-slip boundary condition at the 
fluid/solid interface can be achieved by using a standard 
method of bounce-back on links (BBL) [TH [H]. When 



2 




Pusher 


15(d) 

, / * / \ /-VC f 


'/ ')-25q 




0 5^10 15 20 25 
10^ LB steps 


^-304: 
35^ 
40 


FIG. 1. (a) A cartoon showing an example of a squirmer near 
a flat wall, defining the gap size h and angle </> used in the text. 
In the bulk, the squirmer would move with speed uo along the 
direction shown, at an angle (f> from the horizontal, (b) There 
is a repulsive interaction between the wall and the squirmer. 
Examples of steady state (j>{t) and h{t) observed for (c) puller 
and (d) pusher dynamics near a flat wall in the absence of 
thermal noise are also shown. The interaction range is 5c ~ 
0.16J? (dot-dashed line) and the potential diverges at ~ O.lli? 
(dotted line). 


the boundary is moving {e.g. a colloidal particle) the 
BBL condition must be modified to take into account par¬ 
ticle motion [501 • These local rules can include additional 
terms, such as a surface slip velocity (Eq. [^: in this way 
it is possible to simulate squirming motion mull]. Our 
implementation also includes thermal noise [23] , allowing 
for simulations with a finite Peclet number. 

Simulation parameters: We limit our simulations to 
simple squirmers with = 0, n > 3, but consider both 
pushers (i ?2 < 0) and pullers (i ?2 >0). In simulation 
units (SU) we measure the lengths in lattice spacings 
and time in simulation steps. Parameters, all given in SU 
are: Bi — 0.0015, B 2 = ±0.0075, (which gives a swim¬ 
ming velocity in the bulk equal to mq = = 10“^ 

and j3 = ^ = ±5), fluid viscosity ry = 0.1 and thermal 
noise ksT — 10“®. We considered a fully resolved swim¬ 
mer with radius R = 9.2 (Fig. 1(a)). In order to model 
wall-particle repulsion, we employ a soft potential, V(h), 
which goes smoothly to 0 as the wall-particle separation, 
h (the gap size), approaches (5c, and diverges as h ^ 0 
(see Fig. 1(b), and Supporting Information [24] ). The 
physics is governed by two main hydrodynamic dimen¬ 
sionless quantities: the Reynolds and Peclet numbers. 
Using the parameters above, these are Re=^^^ « 0.09 
and Pe = « 2 x 10^ respectively, where Dr = , 

is the rotational diffusion constant. Our simulations were 
carried out in a rectangular simulation box 120 x 120 x 96, 
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FIG. 2. Simulation results when the range of the soft repul¬ 
sion is varied 5c ~ 0.16i? (dotted line), 5c ~ 0.22R (solid line) 
and 5c ~ 0.27i? (dashed line), with Pe ~ 2 x 10"^. Time devel¬ 
opment of h{t)/R (a) for a pusher and (b) for a puller as well 
as <)>(t) for (c) a pusher and (d) a puller. (The sudden change 
in behaviour for a puller is discussed in the text). The steady 
state h and cf) as a function of 5c for (e) a pusher (the error- 
bars gives the amplitude of the oscillations) and (f) a puller. 
(Initial conditions: ho ~ I.IR and (f)o = 45.0°; Pe « 2 x 10^). 


with periodic boundary conditions in X and Y and solid 
walls at 2 = 0 and z = 95. To see how SU relate to 
physical units, we can, e.g., map a single length and time 
SU to ^ 1/rm and 0.1/xs, respectively [5K] . 

Results: In Fig. [^c,d) we plot the evolution of the 
dimensionless gap size e = h{t)/R and of the angle (j) 
between the squirmer direction and the surface plane 
(Fig. [^a)). For a puller, after an initial collision with 
the soft repulsive wall, the hydrodynamically induced 
torques rotate the particle so that it settles to swim par¬ 
allel to the no-slip wall (beyond the excluded volume 
interaction range; dot-dashed line in Fig. [^c)) with a 
distance h ~ 0.2R, in very good agreement with theo¬ 
retical predictions m- In steady state the puller points 
towards the surface, (() ~ 24 degrees (Fig. [^c)). For a 
P = —5 pusher, previous theories based only on hydro- 
dynamic interactions predict no stable swimming near a 
surface m- However, in experiments, phoretic swim¬ 
mers which are thought to be pushers for mechanistic 
reasons [26], are typically observed to accumulate and 
undergo stable swimming at no-slip surfaces [5113]. Strik¬ 
ingly, our simulations (Fig.[^d)), show a stable periodic 
orbit both in h{t) and 4>(t). During a collision with the 
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soft-repulsive wall, hydrodynamic torques reorient the 
pusher (Fig. [^d)) leading to it swimming away from the 
wall (cj) is on average < 0 in Fig. my This much is 
expected; surprisingly, the long range hydrodynamic in¬ 
teractions between the swimmer and the wall, lead to 
another reorientation of the pusher so that it starts to 
swim towards the wall again. The cycle repeats lead¬ 
ing to the the hydrodynamic oscillations near the no-slip 
surface (Fig. Hd)). While (j){t) < 0 during the oscil¬ 
lations, significant part of the trajectory h{t) is spent 
beyond the external repulsion range (dot-dashed line in 
Fig. m))- Experimentally, these oscillations would be 
difficult to distinguish from true, steady-state trapping, 
providing a potential explanation for the experimental 
observations. Interestingly, for both pusher and puller 
dynamics, the hydrodynamically induced attraction is 
strong enough to resist the effects of thermal noise (see 
Fig. m ,b)). Decreasing /3 reduces the strength of hy¬ 
drodynamic torques |4], we observed no trapping for 
/3 = 0, ±2 (see SI, Fig. SI). 

The external soft repulsion, and in particular its range, 
plays a key role in determining the swimming dynam¬ 
ics. This can be seen from the and h{t) curves 
presented in Fig. [^a-d), for different repulsive ranges 
{6c — 0.16i?, 0.22i? and 0.27i?): these simulations in¬ 
clude the effect of thermal noise, with Pe « 2 x 10^, and 
were all initialised with Hq « I.IR and = 45°. For 
all ranges considered, both the pusher and the puller are 
found to swim near the surface (Fig. [^a,b)). However, 
the hydrodynamic oscillations in the pusher dynamics 
(visible both in h{t) and (j){t) Fig. [^a,c)) are suppressed 
when the repulsive range is increased, and disappear al¬ 
together for 5c = 0.27i? (Fig. [^. In this case, for both 
pusher and puller steady state h < Sc- The plot of the 
swimming orientation ^(<) (Fig. [^c,d)) confirms the ab¬ 
sence of oscillations: the pusher swims by keeping a sta¬ 
ble orientation tilted away from the wall, with (j) slightly 
decreasing when 6c is increased (Fig. [^e)); the puller 
instead is rotated by hydrodynamic torques to point to¬ 
wards the wall, so that (p ~ 90° (this is always the case 
as soon as Jc > 0.22i? (Fig. [^f))). 

Most existing theories of swimmer hydrodynamics rely 
on the far-held approximation which is based on the ve¬ 
locity held a swimmer generates at distances which are 
large with respect to its size. The far-held approxima¬ 
tion can be adapted to include a no-slip wall [T]: as a 
result one obtains the following expressions for the time 
derivative of (/) at a given time. 


d0 Uq 

9,5 sin(20) 


3cos0 f 

1 

64 

Uv 

16 Uv 


where ^ and h' is the distance from the centre of the 
particle to the conhning wall {h' = h + R in Fig. 
Alternatively, one may use lubrication theory to com¬ 
pute the following prediction, based on near-held hydro- 
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FIG. 3. Observed simulation (solid line), far field (dashed 
line) and combined near and far-held (dotted line), results 
for Uci,{t) / {uo /R) for (a) pusher and (b) puller. The repulsive 
range was Sc « 0.16i? (initial conditions hg ~ 3.25R and 
00 = 10°; no thermal noise). 


dynamics [2^ 12^ . 

where e = {h — R)/R [24] . 

To allow detailed comparison between our model and 
the theoretical predictions (Eqs. 2 and 3), we carried out 
simulations starting with hg « 3.25i? and 0o = 10° ^^nd 
at each point we calculated the expression for dp/dt ei¬ 
ther directly from our numerics, or by substituting the 
instantaneous values of h{t) and 0(t) into Eqs. and 
these two equations respectively provide the far- and 
near-held estimate of the system evolution given its cur¬ 
rent state and can be combined by means of a matched 
asymptotic expansion [24] . 

For early times, there is good agreement between the 
rotational dynamics, dp/dt, predicted by the far-held ap¬ 
proximation and that found in our direct numerical sim¬ 
ulations (Fig. [^: we observe a decrease (pusher) and 
increase (puller) of 0(t) from the initial 0o = 10°. Later 
on, the far-held estimate no longer captures the dynam¬ 
ics observed in simulation. In steady state, the far held 
predicts ^ > 0 while simulations show no net motion 
of 0(t), as shown in Fig. [^c,d). When incorporating the 
near-held contribution, we observe very good agreement 
between the theory and our simulations, including the 
trapping of the puller and the oscillations of the pusher 
(Fig.§. This result can be understood by noting that 
the far- and near-held contributions are qualitatively dif¬ 
ferent. In the far-held, a pusher swims stably parallel to 
the wall, whereas a puller rotates until it is perpendicular 
to it [T|. In the near-held, it is the puller which swims 
stably along the wall, pointing slightly towards it m, 
whereas the pusher has no stable swimming solution. 
In our simulations, the particle is trapped close to the 
wall, so near-held hydrodynamics dominates, although 
far-held contributions are non-negligible. [24| An analy¬ 
sis of the dynamics of approach to the surface, dz/dt [53| , 
leads to similar conclusions: prior to interacting with the 
repulsive wall, the far-held works well; when the repulsive 
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FIG. 4. (a) Simulations and far field results for the observed 
steady state swimming speed along the wall U||(t)/uo for a 
pusher (/3 = —5); simulations (solid line) and far field predic¬ 
tions (dot-dashed line) and puller {j3 = -1-5); simulations (dot¬ 
ted line) and far field prediction (dashed line). {Sc ~ 0.16i?; 
initial conditions: ho « I.IR, (j)o = 45°; no thermal noise), 
(b) Time averaged (u| |)/uo as a function of the repulsion range 
Sc, from simulations (open symbols) and far field calculations 
(closed symbols), for both pushers (squares and upward tri¬ 
angles) and pullers (circles and downward triangles). Initial 
conditions: ho ~ 1.17?, tpo = 45°; Pe « 2 x 10^. 


interaction is reached, there is a notable disagreement 
(see SI, Fig. S2). 

For movement along the wall {dx/dt = M||), the sim¬ 
ulations and far field predictions agree reasonably well 
at all times (Fig. Qa); similar conclusions were reached 
by Spagnolie and Lauga who studied the dynamics of a 
swimmer before collisions with the wall my The steady 
state velocities for both swimmers are considerably larger 
than in the bulk, i.e. the presence of a surface accelerates 
the motion (~ 50% increase, see Fig. [^a)). The increase 
in the swimming speed near a solid surface can be un¬ 
derstood intuitively by considering the swimming mech¬ 
anism. The pusher is propelled from behind, thus when 
pointing away from the surface it is pushing the fluid 
flow against a solid wall - this should enhance the swim 
speed, as predicted for swimmers in porous media HES]. 
The speed increase is retained for the periodic swimming. 
The speedup of the puller can be understood in a simi¬ 
lar way: the squirmer is now oriented towards the wall 
so by pulling inward along its swimming axis, it pulls 
itself along the wall. Recent experiments reported an en¬ 
hanced swimming speed up to ^ 2uo for Janus colloids 
on a water-air interface |31| . 

Increasing the range of the repulsive interaction leaves 
the pusher dynamics along the surface mostly unaffected: 
we find that Ux > uq for all the interaction ranges con¬ 
sidered, as shown in Fig. [^a,b). The case of the puller 
is very different, as any repulsive interaction extending 
past the equilibrium swimming distance ~ 0.27? leads 
to hydrodynamic torques orienting the particle towards 
the wall (Fig. [^d)); see also SI, Fig. S3, which shows 
that the far-held approximation for the tangential speed 
remains good in this case as well). The reorientation oc¬ 
curs independently of initial conditions (in Fig. Qb) the 
initial angle is almost tangential, (j)o = 10°), and leads 


to a dramatic slowing down of the particle, whose mo¬ 
tion virtually comes to a standstill when Sc > 0.227? (see 
Fig.gb)). 

In all the cases we have considered, the rotational mo¬ 
tion has a fundamental role in the dynamics of the par¬ 
ticle, and this is affected by Sc- The soft repulsion only 
slows down the particle movement along the surface nor¬ 
mal (as visible from SI, Fig. S2 for dz/dt), and it does 
not create any torques; therefore any rotational motion 
of the particle only arises from the combination of hydro- 
dynamic and Brownian forces. 

Conclusions: We have presented a study of fully re¬ 
solved spherical squirmers swimming between two solid 
walls, using a microscopic model which prescribes a slip 
velocity at the particle surface. Our results show that 
repulsive interactions, which have been neglected in pre¬ 
vious theories of swimmers interacting with surfaces, play 
a very important role in the squirmer’s dynamics. First, 
they can stabilise hydrodynamic oscillations of a pusher 
close to the wall. A recent systematic investigation has 
demonstrated that in the parameter range we consider 
(/3 ^ —5 or below) there is no stable bound state with 
the pusher swimming near the wall m- While experi¬ 
ments routinely observe that bacteria (which are known 
to be pushers) or phoretic swimmers (which are thought 
to be pushers) are attracted to and swim near flat sur¬ 
faces Hi- One way to reconcile these results is if the 
trajectory of the swimmer at late times is oscillatory (a 
limit cycle in the {h, cj)) plane) instead of having con¬ 
stant velocity (a stationary point in the {h,(j)) plane). 
While this conclusion should hold qualitatively for sev¬ 
eral different pusher swimmers, we note that a spherical 
squirmer model does not provide a quantitatively accu¬ 
rate description of a rod-like bacterial swimmer such as 
E. coli, so that the details of its hydrodynamic oscillations 
may in practice differ from those presented here. 

Second, we find that the swim velocity of a pusher is 
much increased with respect to the bulk limit: this be¬ 
haviour can be understood as the particle, on average, is 
directed away from the wall and pushes on it, enhanc¬ 
ing its speed. Third, we find that the tangential velocity 
of a puller slows down dramatically with the range of 
the repulsive interaction with the wall. Our results crit¬ 
ically require near-held hydrodynamics, as the far-held 
approximation poorly captures the rotational dynamics 
we observe. 

Our hndings further imply that the existence and ex¬ 
tent of steric or electrostatic repulsion of the wall could 
be tuned to control properties such as the number density 
and speed of active particles near a surface. Experimen¬ 
tally this could be achieved, by varying either the buffer 
concentration (for electrostatic repulsion) or the polymer 
coverage of the surface (for steric repulsion). These pre¬ 
dictions should be testable with experiments using bacte¬ 
rial swimmers or artihcial microswimmers, although for 
phoretic particles one may need to hrst estimate the ef- 
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feet of chemical gradients, here neglected, on the dynam¬ 
ics IS191132H34]. 
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Supplementary material for hydrodynamic oscillations and variable 
swimming speed in squirmers close to repulsive walls 


SOFT REPULSIVE POTENTIAL AT THE WALL 

We use the following potential between the particle and the wall, 

VUh) = V{h) - V{K) -{h- K) , (1) 

where the wall-particle-surface separation h = ± z — R, and 

Vih)=e{a/hr. ( 2 ) 

Vcs has been cut-and-shifted to ensure that the potential and force go smoothly to zero at h = he- Parameters were 
chosen as e — 0.004, cr — 0.1, zz — 1.0. By choosing ^max — ^top (^c ^c) (^min — ^bottom (^c ^c)) and 

keeping he = 0.5SU constant, we can have a well defined repulsion range 6c, while keeping he and thus the potential 
form constant (Fig.1(b) in the main text). For the calculation of the gap size between the squirmer and the solid 
surface h (Fig. 1(a) in the main text), we define the wall location half-way between the solid node and first fluid node 
(^bottom = 0.5 and Z^op = 94.5), as customary in LB simulations. 


EAR-FIELD APPROXIMATION 


The far-held approximation is based on the velocity held which a swimmer generates at distances which are large 
with respect to its size. The far-held approximation can be adapted to include a no-slip wall [T]: as a result one 
obtains the following expressions for the time derivative of the positions parallel x and perpendicular z as. 


dx 


cost)) 


9Psm{2(j)) f R\^ coscj) f R^ ^ 


32 


h' 


h' 


dz 

dt 


= uo 


smffl — 


9,0 (I - Ssin^ ())) f R\^ sin 


32 


h' 


(3) 


where P = ^ and h' is the distance from the centre of the particle to the conhning wall {h' = h-\- Rva Fig. 1(a) in 
the main text). 
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FIG. SI. Observed simulation results for h{t)/R for a neutral squirmer (/? = 0), for a pusher (/3 = —2) and for a puller 
(/3 = +2). The repulsive range was 5^ ~ 0.16i? and Pe « 2 x 10^. 



3.5 
3 

2.5 
2 

1.5 
1 

0.5 

0 


DC 

■c 



3.5 
3 

2.5 
2 

1.5 
1 

0.5 

0 


DC 


10 LB steps 


10 LB steps 


FIG. S2. Observed simulation results (solid line and dot-dashed line) and far field results (dotted line) for Uz{t)/uo for (a) 
pusher and (b) puller. The disagreement between simulations and far-held approximation persists after correcting the latter to 
include the repulsive interaction from the wall, via an extra normal velocity equal to Uw = — - ^ (where 7 = GirriR). The 
repulsive range was 5c ~ 0.167? (initial conditions ho ~ 3.25i? and (j)o = 10°; no thermal noise). 
















FIG. S3. Simulations and far field results for the velocity along the wall Ux{t)/uo for pusher (/3 = —5) and puller (/3 = +5) 
as well as for a puller (/3 = +5), when the repulsion range was Sc ~ 0.22i? (initial conditions ho ~ S.25R, 4>o ~ 10°; no 
thermal noise). 
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FIG. S4. Definition of the geometry of a spherical squirmer of radius R next to a plane surface. The y direction is into the page. 
In the bulk, the squirmer would move at speed uo along the direction shown, at angle 4> from the horizontal. The squirmer 
is axisymmetric around this axis. Uc is the slip velocity at the point of closest approach, indicated by x. At this point, the 
squirmer wall gap is eR. The fluid flow generates a force and torque Ty on the squirmer, around its centre, which generate 
translation at speed Ux along the x-axis, and rotation, at angular velocity fl around the y-axis. 


LUBRICATION RESULTS 


In Ref. [ 2 ], the wall-parallel force and torque Ty on a squirmer moving near a no-slip boundary is calculated. 
We repeat these calculations in Section]] We obtain different prefactors for and Ty compared to Ref. [2], but we 
agree as to the functional form. In summary, our corrected results for a squirmer oriented at angle (j) away from the 


parallel to the plane (see Fig. S41 are 


F. = -^^l0g(l/e) + O(l), 


Ty = 


5 

IbTTTyR^Uc 


log(l/e)+0(l). 


( 4 ) 


( 5 ) 


Here, Uc is the surface slip velocity (parallel to the wall) at the point of closest approach (x) between the squirmer 
and the wall, which, for the squirmer defined in Eq. I of the main text is 


Uc = —u(7r/2 — (p) = ~ 2 ^o cos (j){l + P sin(()) . 


( 6 ) 


From standard results for the drag on a sphere near a wall [ 3 ] the force and torque give simple expressions for the 
total rotation H and speed of the squirmer 


or, in terms of uq and (p 




n T 

Uy; = — = 0 + 0(1/loge). 




( 7 ) 

( 8 ) 


dx 

dt 


= 0 + 0 


1 


loge 


(9) 


In other words, the term of order unity in the total translational motion of a squirmer near a wall vanishes, and 
the leading order speed decays as 0(1/loge) as the squirmer approaches the wall. It is not possible to calculate the 
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numerical value of this term from lubrication theory, since it depends on longer-range interactions between the whole 
squirmer and the wall [3] . Since this logarithmic decay is very weak, the leading order term will remain comparable to 
uq except for squirmers extremely close to surfaces. In the current simulations e ^ 0.1, giving |l/loge| ^ 0.4, which is 
not small. Hence, it is not contradictory that, in simulations, we see increase as the squirmer approaches the wall: 
this is probably because the swimmer does not approach the wall very closely in the simulations. The lubrication 
calculations merely predict that, for a sufficiently close approach, the translational speed of the squirmer will begin 
to decrease and eventually slow to zero. For the vertical speed, Uz, the term of order unity also vanishes [5], so we do 
not calculate Uz here. For the rotational motion, the next-to-leading order term also decays as 0(1/log e), so it will 
also be significant. 

We can provide an intuitive justification for Eq. 0 -(§. As the swimmer gets closer and closer to the surface, 
most of the viscous dissipation will occur in the thin region around the contact point. We would therefore expect 
the solution to minimise the dissipation in this region. This can be done by ensuring that there is no difference in 
fluid velocity between the particle and the plane surface at this point of contact. Hence, the total velocity on the 
particle surface, taking into account the slip velocity and the solid-body motion of the particle should, in the limit of 
infinitessimal gap size, approach zero, i.e., 

lim {uc + Ux — AH) = 0 . (10) 

€—>■0 

This condition is satisfied (but not uniquely) by Eq. 0-(|8]). It is not satished by the original result derived in Ref. g]. 


CALCULATION OF LUBRICATION FORCE AND TORQUE 


We briefly repeat here the lubrication calculations of Ref. [5] , to obtain the results in Eq. This calculation 

is identical to the standard calculation of the forces and torques of a no-slip sphere near a surface [3], except for 
the new boundary condition on the sphere surface introduced by the finite slip velocity. We hrst define a cylindrical 
coordinate system with = x'^ + y'^ and tan-^ = yjx, where the origin of the coordinate system is the 

point on the plane immediately above the squirmer’s centre. The boundary of the squirmer is defined hy z = h{p*, (/), 
and in the vicinity of the contact point is given by 


h = 


-R 


2A2 


O 



( 11 ) 


To ensure that the equations of motion are all of order unity, we use the dimensionless stretched variables 
A, Y, Z, H, p, with the scaling 


e^/^RX = x, e^/'^RY = y, = p* 

eRZ = z , eRH = h. 


( 12 ) 

(13) 


The stretched height H is 

H = -l-^+0{e). 


(14) 


The fluid velocity field u has cc, p, z components u, v and w respectively. In the stretched coordinate system, the 
Stokes equations are 




\ 


u = A f ei/2_^. 

dX 


l/2_^ ^ 

BY' dZ 


A/2 


du du 
1)X^ dY 


dw 


= 0 . 


(15) 

(16) 


where Vy = xd/dX + ydjdY and with fluid viscosity rj and pressure p. No-slip boundary conditions apply on 
the plane surface: {u,v,w}\z=o = Oj and we write the boundary velocity on the upper surface as u{Z = H) = U, 
v{Z = H) = V and w{Z = H) = W. Expanding the boundary conditions as power series in orders of around 
p = 0, we have 


C/=Uc + 0(ei/"), 
U=0-kO(e^/^), 

W =0 — pcos tjj + 0{e). 


(17) 

(18) 
(19) 
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Performing a similar expansion for the bulk velocity and pressure gives 

u = uo + + 0{e ), 

V = Vo + 


p = e-3/2 


0(e), 

w = 0 + e^^'^wo + ewi + 0{e 
Po + + 0{e) 


3/2 


), 


From Eq. (15), po is independent of Z. Solving for uq and vg in Eq. (15)-(16l then gives 


R dpo 


Z 


Uo = — -;P-{Z - H)Z A - 

2pdX^ ’ iJ ’ 

^ - H)Z. 

“ 2pdY^ ’ 

Combining these solutions and using the equation of continuity (Eq. @) yields, after some algebra 


12p 


V||Po - 


H^R 


4^ p • V||Po - 2“c/^cos'0 = 0. 


Inserting the ansatz po = qoip) cosip then gives an equation in terms of p alone 


H^R d f dqo 


I2p^ri dp 

which has the particular solution 


dp 


H^R 

Vlp^rj 


qo 


H^R dqo , 1 
Ap dp 


-Uc = 0. 
2 


6pp 


As discussed in [5] , the conditions that po be finite everywhere means that this is the only physically relevant 
Next, we rewrite the x,y velocities in the cylindrical polar coordinate system, i.e.. 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


(24) 

(25) 

(26) 

(27) 

(28) 
solution. 


Using the ansatz 


we obtain for the in-plane components 


U = UpP + U^l\! -1- UzZ . 

(29) 

Up = COSIpUpip) , 

(30) 

u^ = sin ipu^ {p), 

(31) 

Up = ^{Z - H)Zq'g + ^Ua, 

(32) 

= -x-(^ - H}Z - —Uc , 

2p pH 

(33) 


where the prime indicates the radial derivative djdp. 

To obtain the total horizontal force on the swimmer, we integrate small elements of force over the swimmer 
surface S', i.e.. 


where [3] 


F,= [ dF, , (34) 

JS 

dFx = X • cr • fids, (35) 


with ff the stress tensor, and dS an infinitesimal area element. We evaluate the stress tensor in cylindrical polar 
coordinates, but use spherical polar coordinates centred on the particle centre, with the polar angle x = 0 at the point 
of closest approach, to specify the normal h. This gives 

dUj, = [—p sin X COS')/'+ (sin X cos— sin X sin + cosy cos — cosy sin V'T.j/iz)] dS. (36) 
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where r is the rate of strain tensor, with components (in the unstretched coordinates 


T -2^ 

op^ 


'T'ljjp^ — p ^ 


d 




1 dVn 


dvo . dv. 


'^P^Z - r, “t“ 


dp* \p* J p* dip ’ 

1 dvz 


dz dp* 


_ dv^ _ 

dz p* dip ’ 


(37) 


and dS' = sin^d^d^/" is the area increment in spherical polar coordinates. Inserting the expansions for the velocities 
and rescaling into the stretched coordinates gives 


dif, = -e sin x cos^ ^ ^ 


2e sinxcos^ — e sinxsin^ ?/'( P 


dp 


d ^ U'tp 

dp 


1 dup 
p dip 


_idup du^ 


dZ dp 
Performing the integral over ip gives 

-^-3/2, 


+ cos X cos^ Ip I e ^ ) — cos x sin^ ip [t ^ - 


du 


ij 


1 du. 


dZ p dip 


dS 


(38) 

(39) 


Ft = R TT 


9o sinx 


-1^ 

dZ 


2e ^/^sinx 

dup 

dp 

e ^/^sinx ( 

'4i 

V dp \ 

'P duz \ 

- cos X 

/ _ 2 du'tp 

1 duz 

? ■ 


p dip 


Z^H 


ti-0 


1 dup 
p dip 

sin xdx 


(40) 


Z^H 


To perform the integral over Xj we expand to first order around x = 0, giving sinx = e^^^p + 0(e). The inner, 
lubrication region extends to some real distance pg of order the particle size, pg = DR, where D = 0(1) is an 
unknown constant which can be obtained by matching to the outer solution. In the stretched coordinate system, the 
corresponding limit is pg = De~^/'^. To lowest order 


Fj; = R^-k 

and evaluating this integral gives 


-9oP 


2 up / du 


R\dZ 


Z=H 


du^tf, 


Z=H 


dp. 


^^. = -^log(po)+0(l). 


(41) 


(42) 


We wish to express in terms of e. As discussed in Ref. HE], we do not need to determine the unknown constant 
D in order to do this, because logpg = —(1/2) loge + logF, so the value of D can be absorbed into the 0(1) term. 
Hence we obtain the expression in Eq. @- 

The torque Ty can be calculated in the same way by integrating the differential elements of torque H 


dTy = i? [(n • z)x — (n • x)z] • cr • ndS”, 
giving, after the same steps as above, the integral 


Ty = R'^tt 


VP ( dup 
R I dZ 


Z=H 


du^ 


Z=H 


dp. 


(43) 


(44) 


which yields the expression in Eq. ([^. 


MATCHING LUBRICATION AND FAR-FIELD RESULTS 


In order to obtain a result which can be compared with the simulation results everywhere, we perform a matched 
asymptotic expansion of the near-held and far-held results. We dehne q = R/h' = 1/(1 -f e). Then, the far-held 
corresponds to g —>■ 0, while the near-held corresponds to e —>■ 0. In the near-held, the next-to-leading-order term is 
0(1/loge), so, in order to match this term to the far-held we dehne the function 


/( 9 ) = 



(45) 
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which has the near-held limit / —>■ —2/log(e), and the far-held limit / —)■ 1. For intermediate values, f{q) is smooth 
and monotonic. We then use the following matched expansion 


1 d(/. 

Ruo dt 


cos (j)q^ 


3 

2 



+ ciq^ 



-|- /3sin {2(j))q^ 


3 

4 


(s + '“"O 


(46) 


where ci and C 2 are constants to be determined by matching to the simulations. This expansion matches both the 
lubrication results and the far-held results in their respective domains of applicability, with corrections of 0(1/log e) 
in the near held, and 0{Pq^) and 0{q^) in the far-held, which is the next order of approximation there [T]. There are 
two next-to-leading-order terms in the far-held because we have linearly independent contributions from the n = 1 
and n = 2 Legendre components of the slip velocity. 

With the htting parameters, ci = 1, C 2 = 0.29, we obtain semi-quantitative agreement with the simulation results. 
Because of the very slow decay of the next-to-leading-order terms in the lubrication theory, we would not expect an 
exact match. Thorough testing of the lubrication theory would require simulations where the swimmer approaches 
much closer to the plane surface. 
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